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We study the chiral properties of domain-wall quarks at high temperatures on an ensemble of quenched con- 
figurations. Low lying eigenmodes of the Dirac operator are calculated and used to check the extent to which the 
Atiyah-Singer index theorem is obeyed on lattices with finite A^s. We calculate the connected and disconnected 
screening propagators for the lowest mass scalar and pseudoscalar mesons in the sectors of difi'erent topologi- 
cal charge and note that they behave as expected. Separating out the would-be zero eigenmodes enables us to 
accurately estimate the disconnected propagators with far less effort than would be needed otherwise. 



1. Introduction 

We use Shamir's formulation [Q of lattice 
domain-wall fermions, which have exact chiral 
flavour symmetry for infinite 5th dimension {N^), 
to study the chiral properties of quarks in high 
temperature QCD and their relation to topology. 

The chirally symmetric Dirac operator should 
have a zero mode associated with each instan- 
ton, and obey the Atiyah-Singer index theorem. 
In addition, these modes yield the disconnected 
contributions which distinguish the vr from the 
rj' screening propagators and the <j{fo) from the 
(5(ao) propagators, and give related contributions 
to the connected propagators. Since we use 
quenched configurations, we cannot determine 
whether the anomalous U{\) axial symmetry is 
restored in the high temperature phase of the 
massless Nf — 2 theory, but we can check re- 
lations required for the U{1) axial symmetry to 
remain broken. 

We measured eigenmodes of the domain-wall 
Dirac operator as a function of A^s on a set of 
quenched configurations in the high temperature 
phase. We check the Atiyah-Singer theorem for 
Nr, = 10 and calculate connected and discon- 
nected scalar and pseudoscalar screening propa- 
gators. Early results were presented at Lattice'98 
|2|. Section 2 describes our analysis and tests of 
the index theorem. The meson screening propa- 
gators are discussed in section 3. Section 4 gives 
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our discussions and conclusions. 



2. Domain-wall quarks at high tempera- 
tures 

For our study of domain-wall quarks at high 
temperatures we use 16"^ x 8 quenched configura- 
tions 8X13 = 6.2 (170), (i = 6.1 (100) and /3 = 6.0 
(100). (At Nt = 8, /3c « 6.0.) On each configura- 
tion, we estimated the topological charge by the 
cooling method. 

We studied eigenvalue trajectories of the her- 
mitian Wilson Dirac operator 75 -Dm^ as the bare 
mass was varied. These vanish at would-be zero 
modes associated with instantons Based on 
these studies we chose the Shamir domain-wall 
mass parameter M = 1.7 for all 3 (3 values. 

We calculated the lowest 2 eigenmodes of the 
hermitian domain-wall Dirac operator j^TZDdw 
for all f3 = 6.2 configurations with non-trivial 
topology and for all (3 = 6.1 and /3 = 6.0 configu- 
rations, for N5 = 4, 6, 8 and 10. (A third N5 = 10 
eigenmode was calculated for each P = 6.0 con- 
figuration with topological charge ±3.) 

For /3 = 6.2, there is a clear separation between 
would-be zero modes whose eigenvalues decrease 
exponentially with Nc, and non-zero modes which 
rapidly approach a constant. By f3 = 6.0 the 
situation is less clear. For f3 = 6.2, A^5 = 10 
appears adequate. 

We calculate the chiral condensates (qq) and 
(qIsq) for case A^5 = 10, using an eigenmode 
enhanced stochastic estimator. For each configu- 
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ration, exact chiral symmetry requires (9759} to 
obey the Atiyah-Singer index theorem: 



''^{q{x)i5q{x))u = Qu- 



(1) 



As figure |^ shows, the index theorem a.t f3 — 6.2 
is well approximated down to masses comparable 
to the eigenvalue of ■'jsTZDdw at m = of smallest 
magnitude. 



1.50 



1.25 - 



1.00 



0.75 



0.50 



0.25 



0.00 



13 = 6.2 



X — direct calculation 



using Ward identity 




0.00 0.02 0.04 0.06 0.08 0, 
m 



10 



Figure 1. mTr(75 5'_F) as a function of mass for 
(3 — 6.2. The solid lines for Q = 1 are the contri- 
butions from the lowest eigenmode. 



Q = and on each z slice for Q / 0. The discon- 
nected parts of these propagators are measured 
separating out the contributions of the eigen- 
modes (projected on to the domain walls) and ap- 
proximating the remainder of the required traces 
with a stochastic estimator. For a discussion of 
stochastic estimators for staggered quarks seeQ. 

The P — 6.2 disconnected propagators behave 
as expected. For Q = 1, the contribution is size- 
able and increases as quark mass decreases (it 
should diverge as if chiral symmetry were 

exact), as shown in figure 0. For Q = the con- 
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3. Meson screening propagators. 

Screening propagators, i.e. propagators for 
spatial separations, describe the excitations of 
hadronic matter and the quark-gluon plasma at 
finite temperature. We measure the connected 
and disconnected parts of these propagators for 
scalar and pseudoscalar mesons with zero mo- 
menta transverse to their separation (Z). 

The connected parts of these propagators are 
measured with a noisy point source on 1 z slice for 



Figure 2. The disconnected contribution to the 77' 
propagator at (3 — 6.2. The corresponding graph 
for the a is virtually identical. 



tribution is extremely small compared with the 
connected contribution (on the scale of figure ^ 
all points would lie on the horizontal axis). This 
is consistent with the expectation that the Q = 
contribution vanishes in the chiral limit. 

For the connected propagators, we consider the 
combinations ^{P-k ± Ps) (Note: the connected 
part of the rj' propagator is the tt propagator and 
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the connected part of the a propagator is the 5 
propagator). For Q 0, flavour chiral symmetry 
restoration and the absence of a disconnected con- 
tribution would make \ {P'k-\-Ps) vanish in the chi- 
ral limit. The observed value is very small. The 
other combination, \{P-n — Ps), is non- vanishing 
and has a finite chiral limit, as expected. 

In the Q = 1 sector, ^{Ptt + Ps) should equal 
the disconnected part of the a and rj' propagators 
in the chiral limit, for chiral flavour symmetry 
restoration with broken C/(l) axial symmetry to 
be possible for A^/ = 2 |||. i(P^ + Pg) is plotted 
in figure H. Figures I and I are consistent with 
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tial in N^. The index theorem is well approxi- 
mated and the tt, cr, rj' and S propagators show 
the correct behaviour for the restoration of chi- 
ral SU{2) X SU{2) fiavour symmetry for Nf = 2. 
The disconnected parts of the a and rj propaga- 
tors come entirely from the Q sector and have 
the behaviour needed to give a finite contribution 
when the Nf = 2 determinant is included. Isolat- 
ing the lowest eigenmodes greatly improves our 
stochastic estimation of traces. 

Our preliminary results for /? = 6.0 indi- 
cate a more complex behaviour. It is not yet 
clear whether this indicates that — 10 is 
too small, or that the condensation of instanton- 
antiinstanton pairs, as observed by Heller et al. 
for Nf = A using Ginsparg- Wilson fermions |^ , is 
playing an important role. 

A dynamical simulation is needed to determine 
if, for Nf = 2, the U{1) axial symmetry remains 
broken above the transition. Such simulations 
are being performed by the Columbia group at 
= 4 [0, and preliminary indications are that 
it does remain broken. These pioneering results 
are limited by the strong couplings dictated by 
Nf = 4:, and use of heavier quarks than would be 
desirable. To push closer to the chiral limit prob- 
ably requires isolating the lowest eigenmode(s) in 
the simulations. 

These calculations were performed on the C90 
and J90's at NERSC. 
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Figure 3. The connected propagator ^{P-k + Ps) 
for Q = 1 at /3 = 6.2. 



this expectation. 

4. Discussion and conclusions. 

For quenched lattice QCD at A^t = 8, domain- 
wall quarks with N^ > 10 are a good approxima- 
tion to chiral quarks at /? = 6.2 and suggest an 
approach to chiral symmetry which is cxponcn- 
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